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C^ . Abstract. Many astrophysical flows occur in inhomogeneous media. 

CN • The broad-line regions (BLR) of active galactic nuclei (AGNs) are one 

of the important examples where emission-line clouds interact with the 
outflow. 

We present results of a numerical study of the interaction of a steady, 

planar shock / supersonic postshock flow with a system of embedded 

cylindrical clouds in a two-dimensional geometry. Detailed analysis shows 

that the interaction of embedded inhomogeneities with the shock / post- 

^J ■ shock wind depends primarily on the thickness of the cloud layer and the 

q | arrangement of the clouds in the layer, as opposed to the total cloud mass 

t — ■ ■ and the total number of individual clouds. This allows us to define two 

classes of cloud distributions: thin and thick layers. We define the critical 
cloud separation along the direction of the flow and perpendicular to it. 
This definition allows us to distinguish between the interacting and non- 
interacting regimes of cloud evolution. Finally we discuss mass-loading 
in such systems. 
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1. INTRODUCTION 



^ Mass flows are important in many astrophysical systems from stars to the most 

distant active galaxies. Virtually all mass flow studies focus on homogeneous 
media. However, the typical astrophysical medium is inhomogeneous with the 
inhomogeneities (clumps) arising due to initial fluctuations of mass distribution, 
the action of instabilities, variations in the flow source, etc. The presence of 
inhomogeneities can introduce not only quantitative but also qualitative changes 
to the overall dynamics of the flow. 

Active galactic nuclei represent one astrophysical site where a shock / wind 
interaction with a system of inhomogeneities may take place. Practically all 
current models of AGNs agree that the emission-line clouds in BLRs are essential 
for explaining the observed properties of AGNs (Urry & Padovani 1995; Elvis 
2000). However, despite the fact that the nature of BELR clouds is an open 
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question, any self-consistent model of the AGNs should properly account for 
properties of BELR cloud interaction with mass outflow. In particular, such 
self-consistency is important in the context of cloud survival time and cloud 
displacement prior to destruction. Even when clouds are magnetically confined 
and stabilized against evaporation (Rees 1987), disruptive action of outflows 
may prevent cloud survival over the dynamically significant AGN timescales. 

(Klein, McKee, Colella 1994) (hereafter KMC) addressed the problem of 
shock-cloud interaction, providing a detailed description of the dynamics of 
a single, dense, unmagnetized cloud interacting with a strong, steady, planar 
shock. That work gives an excellent introduction to the subject, in particular 
the description of the astrophysical significance of the problem of a shock wave 
interacting with a dense cloud (see also (Gregori et al. 2000)). In this work we 
investigate the general properties of strong shock / supersonic wind interaction 
with a system of embedded clouds and determine the key quantities governing 
the evolution of such systems. 

2. RESULTS 

We have numerically investigated the interaction of a strong, planar shock wave 
with a system of dense inhomogeneities, embedded in a more tenuous and cold 
ambient medium. Our code is the AMRCLAW package, which implements an 
adaptive mesh refinement algorithm for the equations of gas dynamics (Berger &: 
Oliger 1984; Berger & Jameson 1985; Berger &: Colella 1989; Berger & LeVeque 
1998) in two dimensions. We have assumed constant conditions in the global 
postshock flow constraining the maximum size of the clouds only by the condition 
of the shock front planarity. Our results are applicable to strong global shocks 
with Mach numbers 3 ^ Ms <^ 1000. The range of the applicable cloud - 
unshocked ambient medium density contrast values is 10 — 1000. Figure 1 shows 
a case of a strong shock and a supersonic post-shock flow interacting with an 
inhomogeneous system of fourteen identical clouds in regular distribution. 

Cloud evolution due to the interaction with the global shock and the post- 
shock flow has four major phases, namely the initial compression phase, the 
re-expansion phase, the destruction phase, and finally the mixing phase. Each 
image in Figure 1 roughly illustrates each of those phases. 

The timescale we use to define time intervals in our numerical experiments 
is the time it takes for the incident shock wave to sweep across an individual 
cloud. This is called the shock- crossing time, t$c = (2a max / 'v s) , where a max is 
the maximum cloud radius in the distribution. 

We define the cloud destruction time tcD as the time when the largest 
cloud fragment contains less than 50% of the initial cloud mass. Typically, in 
our simulations tcD ~ 2Atsc- 

A simple model for the cloud acceleration during the first three phases, i.e. 
prior to its destruction, can be developed. We find the cloud velocity 



v c (t) = < 



vps[\-(j — ai + a 2 ) J, t < \2t S c 



vps\l-{[- 12j b 1 + —a 1 + a 2 ) ), t<t CD 



(1) 
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Figure 1. Run MIA. Time evolution of a system, containing fourteen 
identical clouds in a regular distribution and interacting with a M$ = 
10 shock wave. Shown are synthetic Schlieren images of the system at 
times 22 t S c, 35 t S c, 50 t S c, 69 tsc- 
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Here, vps is the unperturbed postshock flow velocity. For the cases of infinitely 
strong shocks, i.e. Ms — ► oo, the above coefficients have the values 

ai = 1.83 • 1(T 3 ; a 2 = 1.09; h = 8.51 • 1(T 5 . (2) 

Equation (1) shows that the maximum cloud velocity is not more than 10% of 
the global shock velocity and not more than 13% of the postshock flow velocity. 
Cloud displacement prior to its destruction can be described as 



Lc{t) 



1 , ft (Ol\ 



aoci Inl— — +1 , 0<t<12t SC 

a\ ^t S c y a 2 ' 1) 

_t 12 \ { ' 

r/ rj I c 2 tan~ 1 (-^- ) - c 3 , 12tsc <t<t C D 

\t S c V ^W J 



In the limiting case Ms — > oo the values of the coefficients a± and a 2 are defined 

in (2) and ci = 1.5; c 2 = 103.22; c 3 =10.9; c 4 = 9.43 • 10~ 2 ; c 5 = 112.56. 

We find the maximum cloud displacement prior to its destruction, Lcd, is 

L C d = L c (t CD ) < 3.5a max . (4) 

The results of our model are in excellent agreement with the numerical experi- 
ments. The difference in the values of cloud velocity and displacement between 
analytical and numerical results is <J 10%. 

The principal conclusion of the present work is that the set A of all possi- 
ble cloud distributions can be subdivided into two large subsets A/, thin-layer 
systems, and Am, thick-layer systems, defined as 



A/ : (AxN)max < LcD, 

Am : {AxN)max > Lcd, 



(5) 



where (Axi\f) max is the maximum cloud separation in the system along the 
direction of the flow, or the cloud layer thickness. 

Distributions from each subset exhibit striking similarity in behaviour (e.g. 
Figure 2). The systems containing from one to five clouds, arranged in a single 
layer, exhibit exactly the same rate of momentum transfer from the global flow 
to the clouds. The two fourteen cloud runs M14 and Mu r have different cloud 
distributions, different total cloud mass, different cloud sizes. Nevertheless, the 
rate of the kinetic energy fraction increase during compression and re-expansion 
is different from the single layer cases but is still the same for both fourteen 
cloud runs. Other global properties exhibit the same behaviour. We conclude 
that the evolution of a system of shocked clouds depends primarily on the total 
thickness of the cloud layer and the cloud distribution in it, as opposed to the 
total number of clouds or the total cloud mass present in the system. 

The key parameters determining the type of evolution, are the cloud de- 
struction length Lcd, defined above in (3), and the critical cloud separation 
transverse to the flow d cr it- The latter is defined by the condition that the time 
for adjacent clouds to expand laterally and merge into a single coherent structure 
is equal to the cloud destruction time tcD- It can be described by the expression 
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Figure 2. Time evolution of the global average of the kinetic energy 
fraction (??fcj n )2D for the following runs: 

Ml, M2, M3, A5 - systems with 1, 2, 3, 5 identical clouds correspond- 
ingly, arranged in a single row with constant cloud center separation of 
4.0ao; M14 - system with 14 identical clouds with cloud center sepa- 
ration in a row equal to 4.0ao, separation between rows equal to 7.0ao; 
M14 r - system with 14 clouds of random size and placed in random 
positions. 



We should be able to determine, either from observations or from theoret- 
ical analysis, the thickness of the cloud layer (AxN) max . This determines the 
class of the given cloud distribution, A/ or Am- For distributions from the set 
A/ with average cloud separation (Ayjv) > d cr it evolution of the clouds during 
the compression, re-expansion, and destruction phases will proceed in the non- 
interacting regime and the formalism for a single cloud interaction with a shock 
wave (e.g. KMC) can be used to describe the system. On the other hand, if the 
cloud separation is less than critical, the clouds in the layer will merge into a 
single structure before their destruction is completed, and though the compres- 
sion phase still can be considered independently for each cloud, evolution during 
the re-expansion and destruction phases proceeds in the interacting regime. 

When the distribution belongs to the subset Am it is necessary to determine 
the average cloud separation projected onto the direction of the flow (Axat) and 
compare it against Lqd'- if (Axjv) > Lqd evolution of the system can be roughly 
approximated as of a set of distributions from the subset A/ and the above "thin- 
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layer case" analysis applies. If, on the other hand, (Axn) < Lqd (especially if 
(Ay^v) < d cr it) the system evolution is dominated by cloud interactions and a 
thin layer formalism is inapplicable. 

Finally we consider mass-loading. Here our principal conclusion is that 
mass-loading is not significant in the cases of strong shocks and supersonic winds 
interacting with inhomogeneities whose density contrast is in the range 10—1000. 
In part this is due to short survival times of clouds as well as the very low 
mass loss rates of the clouds even during the times prior to their destruction. 
Therefore, mass-loading in such systems is not likely to have any appreciable 
effect on the overall dynamics of the global flow. 

The major limitation of our current work is the purely hydrodynamic nature 
of our analysis that does not include any consideration of the magnetic fields. 
As was discussed by KMC, cold dense inhomogeneities (clouds) embedded in 
more tenuous hotter medium are inherently unstable against the dissipative ac- 
tion of diffusion and thermal conduction. Although weak magnetic fields, that 
are dynamically insignificant up to the moment of cloud destruction, can in- 
hibit thermal conduction and diffusion and, therefore, stabilize the system of 
clouds, those magnetic fields may become dynamically important due to tur- 
bulent amplification during the mixing phase (see (Gregori et al. 2000) for a 
three-dimensional study of the wind interaction with a single magnetized cloud) . 
We intend to provide a fully magnetohydrodynamic description of the interaction 
of a strong shock with a system of clouds in future work. 

Acknowledgments. This work was supported in part by the NSF AST- 
0978765, NASA NAG5-8428, and DOE DE-FG02-00ER54600 grants. 

The most recent results and animations of the numerical experiments, described 
above and not mentioned in the current paper, can be found at 
www.pas.rochester.edu/~wma. 
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